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Abstract 

The minimal supersymmetric standard model leads to precise predictions of the properties of 
the light Higgs boson degrees of freedom that depend on only a few relevant supersymmetry¬ 
breaking parameters. In particular, there is an upper bound on the mass of the lightest neutral 
Higgs boson, which for a supersymmetric spectrum of the order of a TeV is barely above the one 
of the Higgs resonance recently observed at the LHC. This bound can be raised by considering 
a heavier supersymmetric spectrum, relaxing the tension between theory and experiment. In a 
previous article, we studied the predictions for the lightest GP-even Higgs mass for large values 
of the scalar-top and heavy Higgs boson masses. In this article we perform a similar analysis, 
considering also the case of a GP-odd Higgs boson mass of the order of the weak scale. We 
perform the calculation using effective theory techniques, considering a two-Higgs doublet model 
and a Standard Model-like theory and resumming the large logarithmic corrections that appear 
at scales above and below mA, respectively. We calculate the mass and couplings of the lightest 
GP-even Higgs boson and compare our results with the ones obtained by other methods. 
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I. INTRODUCTION 


Since the discovery of the Higgs boson at the LHC in 2012 DEI. both the ATLAS and 
CMS experiments have made increasingly precise measurements of its mass Mh, mainly in 
the h —)■ ZZ, 77 decay channels. Using ~ 5 fb“^ of data at = 7 TeV and ~ 20 fb“^ of 
data at a/s = 8 TeV, the ATLAS and CMS experiments have measured [31 ID 

ATLAS: Mh = 125.36 ±0.37 ±0.18 GeV, (1) 

CMS: Mh = 125.02 QeV, (2) 

where the quoted uncertainties are statistical and systematic, respectively. The result from a 
recent combination of the measurements from ATLAS and CMS is Mh = 125.09±0.21±0.11 
GeV [5]. 

Low-energy supersymmetry (SUSY) is a highly predictive framework that can accom¬ 
modate the observed Higgs mass and Standard Model (SM)-hke properties in a variety of 
models [6]. These models contain at least an extra Higgs doublet and the observed Higgs 
boson is usually identified with the lightest GP-even state h, with properties that deviate 
from the SM one depending on the mixing with the other neutral scalar states in the theory. 
In the Minimal Supersymmetric Standard Model (MSSM), the Higgs sector reduces at tree 
level to a type H two-Higgs doublet model (THDM), with the mass of the lightest GP-even 
Higgs boson bounded to be smaller than the neutral gauge boson mass Mz- 

This tree-level result, however, is modified by SUSY-breaking effects, receiving large 
radiative corrections from heavy top squarks. In the case of heavy supersymmetric particles 
and nonstandard Higgs bosons, the Higgs boson mass may be determined as a function of the 
stop masses and their mixings, depending only weakly on other SUSY-breaking parameters. 
Models with heavy supersymmetric particles are motivated by the absence of any signihcant 
deviations of flavor or precision measurement observables with respect to the SM predictions. 
Hence, a precise computation of the Higgs mass as a function of the stop mass parameters 
is of signihcant interest. 

There has been much activity in the computation of the Higgs mass in the MSSM in the 
past. The Higgs mass has been calculated by performing hxed-order perturbative calcula¬ 
tions in the MSSM, as well as in effective theory analyses, in which the dominant logarithmic 
dependence has been resummed by renormalization group (RG) methods. For supersym- 
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metric particle masses of the order of the weak scale, an accurate prediction of the Higgs 
mass may be obtained by computing the radiative effects diagrammatically up to a hxed 
order in perturbation theory [M2]. Alternatively, the dominant radiative corrections at a 
given order in perturbation theory may be obtained from effective potential methods, using 
derivatives of the effective potential V{Hi, H 2 ), for values of the Higgs held equal to their 
vacuum expectation values {Hi) = vi, (i/ 2 ) = V 2 [T3] - [T6] . These hxed-order calculations 
have been now carried out up to partial three-loop order [T71 - I20] . 

On the other hand, for heavy supersymmetric particles, the effective held theory approach 
may be implemented by integrating out MSSM particles, considering the induced thresholds 
to the relevant couplings and running them down to the electroweak scale, evaluating the 
ehective potential approximation of the Higgs mass, and, after appropriate corrections, the 
pole mass [2T1 - I25] . It is clear that for low values of the supersymmetric particle masses, where 
the logarithmic corrections are similar in size to the nonlogarithmic ones, the hxed-order 
calculations are expected to lead to the most accurate values. For very heavy supersymmetric 
particles, the logarithmic corrections become very large, the hxed-order perturbation theory 
breaks down, and the RG approach leads to an appropriate resummation of the leading 
logarithmic corrections. In this case, the ehective held theory methods may lead to a more 
accurate determination of the Higgs mass. 

In a previous work [26|, we used ehective held theory (EFT) calculations to compute 
the mass of the lightest GP-even Higgs boson in the MSSM, in the case of heavy stops 
and nonstandard Higgs bosons. A similar approach was also taken recently in Refs. [271- 
[29] . We studied the cases of light and heavy charginos and neutralinos, which can lead 
to relevant radiative corrections to the lightest GP-even Higgs mass. I Furthermore, we 
provided an analytical approximation for the relevant three- and four-loop corrections to 
the Higgs mass that revealed a large cancellation between the dominant and subdominant 
leading-log contributions, leading to a large diherence between our computations and the 
previous partial three-loop calculations discussed above. 

In this article, we perform a similar study for the case of a small GP-odd Higgs mass, 
characterizing a light nonstandard Higgs boson spectrum. In this case, the theory below 
the stop mass scale is a THDM, with the possibility of additional charginos and neutralinos, 
depending on the choice of the gaugino and Higgsino mass parameters; see Fig. This 
approach was hrst detailed in Ref. [30]. The presence of two GP-even Higgs bosons at low 
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FIG. 1. Examples of hierarchies of scales examined in Ref. |26j (left), and in this paper (middle, 
right). 


energies leads to mixing effects in the GP-even Higgs sector that affect the predicted lightest 
GP-even Higgs boson mass and couplings, and therefore modify the Higgs physics at high 
energy colliders. 

This article is organized as follows: in Sec. |^we describe the properties of the low energy 


effective theory, the THDM. In Sec. HI we describe the constraints on this generic framework 


when we assume the presence of a softly broken supersymmetric theory. In Sec. |IV| we study 

we 


the numerical predictions for the Higgs boson masses and mixing angles. In Sec. VI 


describe the approach to the alignment limit and the comparison with the values predicted 


in the hMSSM approach. We reserve Sec. VH for our conclusions. 


II. TWO-HIGGS DOUBLET MODEL 


The most general scalar potential with two complex SU{2)l doublet Higgs helds d*!, $ 2 , 
each carrying hypercharge Y = 1, is 


V = - (m?2<hl<l>2 + h.c.) + y(<hl$i)2 + 

+ A3(<I>l$l)(<I>^<h2) + A4(<I>l<h2)(<I>^$l) 

+ \ y (<I>1<I>2)' + [A6(<I>I<hi) + A7(<I>^<h2)] <I>1<I>2 + h.C. 


(3) 


We assume GP-conservation and for simplicity, will take the coefficients and 

Ay to be real. At the minimum of the scalar potential, the Higgs helds acquire vacuum 
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expectation values 


and we can parametrize them by writing 



( 4 ) 






(5) 


where 0^ is complex and 0^, are real. We choose the Vi to be real and non-negative, with 
the usual relations 


V = \Jvl + V 2 — 246 GeV, tjj = ia.\i (3 = V 2 /V 1 . 


( 6 ) 


After electroweak symmetry breaking (EWSB), Eve physical Higgs bosons remain in the 
spectrum: two CP-even, charged two CP-even, neutral h, P; and one CP-odd, neutral 
A. Minimizing the scalar potential, we can eliminate ml,m 2 , and we have the following 
expressions for the squared masses of A and P^: 


= ^12 - 2^ ( 2-^5 + + \7t13) , 

m^± =m\ + - A 4 ). 


(7) 

( 8 ) 


The squared mass matrix for the CP-even, neutral Higgs bosons in the {$ 1 , $ 2 } basis is 


= 


Ml Ml 


•'''‘12 \ ^2 1^/3 


= 


+ c 


/ll /12 


(9) 


Ain Al'22/ V“W <=s / \/i2 

where = sin/3,c^ = cos/3. Throughout this paper, we will employ similar shorthand 
S 0 = sin 6 ', C 0 = cos 9 ,t 0 = tan 6 * for a generic angle 9. The fij are 

/ll = AiC^ -|- 2 \cpSii + AsS^ , 

/12 = (A3 -|- A4)C/3S^ -|- AgC^ + M-lp ■, 
f 22 = A2S^ -f 2\iCi3Si3 + AsC^ . 


( 10 ) 

( 11 ) 

( 12 ) 


Diagonalizing this matrix, the masses of the physical CP-even, neutral Higgs bosons are 
given by 

1 . .0 . ■ ..... - 1 


ml h = 2 (^^M^ ± a/(T rAd^)^ — 4 det M"^^ = - (^Mh + Ml ± 6 m^^ , 


(13) 
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where the mixing angle a for the neutral CP-even states is 


Cry - 


hm2 + - Ml 


'■22 


25'm? 


Srv — 


V 2 Mf 2 



■\/Sm^(dm,^ + Mfi — M22) ' 

$ 1 ' 



( 14 ) 


and the mixing angle is defined in the range —7r/2 < a < 0. 

We can also rotate to the Higgs basis {Hi,H 2 } [31], 

'<hi\ 

J ■ 

where only one of the scalars receives a vacuum expectation value, (Pi) = v/\/2. In the 
Higgs basis, the CP-even mass matrix takes a similar form, 

^0 o' 




Mir = m, 


+ V 


0 1 


2 ( 9ii 9 i2 
.912 922 


where 


9 11 = + '^2'S^ + 2(A3 + A4 + A5)s^c| -f dAgC^s^ + 4 A 7 s|c^ , 

912 = CgSg{\2s‘jj — AiC^ -1- (A3 + A4 + A5)c2/3) + 3(A7 — + AgC^ — XjS^ , 

922 = (Al + A2 )c|s^ — 2(A3 + A4)s^C^ + A5(s^ + C^) + (A7 — A6)s2/3C2/3 , 

and it follows that in this basis the mixing angle is /3 — a, namely 


(16) 

(17) 

(18) 
(19) 




( 20 ) 


C-jB—a. ^fB—a 
^fB—ot C^—Q 

When the mixing Cfs-a is small, this basis is convenient since the lightest CP-even Higgs 
tree-level couplings are identified with the SM Higgs ones. More generally, the Lagrangian 
describing the coupling of the Higgs bosons to the top and bottom quarks at scales below 
Ms may be parametrized in the following way: 

C = {h, + 6hh) + ep {ht + Sht) tnQm + + e,, + H.c. (21) 


From here it follows that the bottom and quark running masses are given by 

h ) ’ 


hbV C , 5hb , 


x/2 

htv 


mt = —^sg 1 + ^ + 


x/2 


5hf Ahf 


ht 


httg 


( 22 ) 

(23) 
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with the relevant couplings evaluated at the weak scale. Observe that while the corrections 
to the bottom coupling are loop suppressed, they are enhanced at moderate or large values 
of tp and therefore they may take values of order 1 in this regime. On the contrary, the 
corrections to the top coupling are suppressed by both loop and tp factors and therefore 
tend to be small. 

At tree level, the MSSM Yukawa couplings are related to the SM Yukawa couplings by 


Vt yh hi)Cpy y^- h^-cy. 


(24) 


From Eqs. (22)-(23), it follows that these couplings are modihed at one-loop order at Mg 


in the following forms 


ht = 

h = 

hr- = 


yt 


1 


S/3 1 + + A/ 

Vb 1 

C/3 1 + + Afe 

Vt 1 


(25) 

(26) 
(27) 


C/3 1 -h (5 t- + At 

where 5/ = 5hi/hi are the terms without factors of tjs, and A* = {Ahtt^^)/ht [A/, = 
{Ahbtb)/hb, Ar = {Ahrtb)/hr] are ty suppressed [enhanced] terms: 


5 8 2 

^ /(m^^, p) - -gYMiAtI{mi^,mi^,Mi ), (28) 


tyAt 


K 


= I im, , rnt2 , "iff) + hlyAb I (m^^, , p) 


- M2, y) + Sfe/(m^^,M 2 ,/i)J + M2, y) + M2, y) 

+ ^hyMiyl ‘^I{mi^,rrn^,Mi) + ^ cll{mi^, M^, y) + sll{mi^,Mi, y) 


- 2 


("^ti ,Mi,y) + ,Mi,y) 


(29) 


^ - hly^ y) + ^glMiAbI{mi^,mi^,Mi) , (30) 

A 8 

—f ^I (rnq, , y) 

Hi Z ^ o 

-glM2yiy cll{mi^,M2,y) + sll{mi^,M2,y) cll{mi^, M2, y) + M2, y) 

- \9lMiy\^I{mi^,mi^,Mi) + ^ c^J(m^^, Mi, p) + sg/(m^^. Mi,p) 

2 T\/T I ^2 


+ 


Sfc ^ ("ifoi, ^ 1 ,/^) + Cft / (m^^, Ml, p) 


(31) 
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( 32 ) 



(33) 


In the above expressions, k = (l/ 167 r^) is a loop factor; At (^fe,r) are the trilinear conplings 
of the stops to the Higgs held $2 (*^ 1 ); are the sfermion eigenstate masses; Mi 2 are the 
hypercharge and weak gangino masses; rrig is the glnino mass; and fi is the Higgsino mass 
parameter. The parameters St,Sh,Sr {ct,Ch,Cr) are the sines (cosines) of the stop, sbottom, 
and stau mixing angles, and the fnnction I[a, b, c) is dehned as 


a?h'^ log(a^/ 6 ^) + b‘^c^ log( 6 ^/c^) + \og{c^/a?) 
(a^ — b‘^){b‘^ — c^)(a^ — c^) 


/(a, b, c) 


(34) 


We will assume that the masses nig = my = my = my = my = Ms (such that = 
= 1/2 with X = t,b,T). We will consider the two scenarios M 2 = Mi = /x = Ms 
and M 2 = Ml = /x = 200 GeV. With these choices, the above expressions contain the 
dominant contributions to the threshold corrections, which also include all terms necessary 
for consistency with our threshold corrections to the quartic couplings. [27] 

Strictly speaking, below M 5 , the couplings Xht^b and ht^ + evolve in slightly differ¬ 
ent ways. However, since the dominant contribution from QCD in the RG evolution of the 
couplings is the same, and the couplings Xht^ are already loop suppressed, we shall approx¬ 
imate the ratios ^ as constants below Ms and concentrate only on the RG evolution of 
the top- and bottom-quark couplings to the helds Hu and respectively. We expect this 
approximation to have a negligible impact on the Higgs boson masses. 


Using the above expressions, one can easily prove that the couplings of the light phys¬ 
ical Higgs boson h to top and bottom quarks and vector gauge bosons are given by (see, 
e.g. Ref. [M] i 



ghVV — Sg-a ghVV y 
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( 35 ) 



















where ghtt^dh^ 9h^v denote the SM couplings of top quarks, bottom quarks and weak 
gauge bosons to the Higgs. One observes that in the regime of moderate or large values of 
tjs and small values of C/ 3 _q,, the bottom coupling can get sizable departures from the SM 
value, while the top and vector gauge boson couplings tend to be close to their SM values. 

Similarly, for the heavy Higgs boson H, one obtains 


gmt — 


gmb — 


Sp-c 


+ CjS-a I + 


('5/3—0 tf) "t" ('13—a) 


f t/S Sj^—c 

1 + (5t + 

Afe ( ^(3 ^13—c 


1 + 5b + Afe 




9htt, 


SM 
ghbb ) 


(36) 


gHVV — Cp-a 9hVV 


Hence, one observes that for small values of c^_q, the coupling gnbb of the heavy CP-even 
Higgs to bottom quarks is affected by loop corrections and can become sizable at large values 
of tp. The top-quark coupling to the heavy Higgs instead remains suppressed by either loop 
or tp factors. 

For completeness, we stress that there is a close connection between the coupling of the 
heavy CP-even Higgs and of the CP-odd Higgs to top and bottom quarks. These CP-odd 
Higgs boson couplings are given by 


gAtt — 


1 


tp Ai 

(1 + 5i + Aj)s^ 


gm 


gAbb = 


t /3 Ab 

(l + 5b + Ab)s2 


9h^ 


(37) 


III. THE MSSM HIGGS SECTOR 

The MSSM Higgs potential for the two Higgs doublets Hd^Hjj with respective hyper¬ 
charges Y = — 1,1 is 

Vh = \{9l + gl)(\HD? - \Hu\^f + \gl\HlHu\^ + \p.\'‘{\Ho? + \Hu?) 

+ + m 22 \Hu? + ■ Hu + h.c.), (38) 

where Ho ■ Hu = eabH^Hu. These originate from the D-terms in the superpotential and 
the soft supersymmetry-breaking terms. To recover the form of the THDM potential in Eq. 
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Q, let ml = mli^ + for k G {1,2} and m^2 = B/j,, with the following relations between 
the helds, 

$1 = —ia2H*j ^, $2 = Hjj. ( 39 ) 


The terms in Eq. (38) become 


^ |<l>in<l>2p - (<hl$2)(<h^^>i), Hd-Hu^ -$I<h2, (40) 


and we have the following tree-level relations for the qnartic couplings: 


Al — A 2 — -{9I - I - 9y) , 

(41) 

-^3 = -^{92 ~ 9y) 1 

(42) 

X 1 2 

-^4 — 2 ^^ ’ 

(43) 

-^5 = -^6 = -^7 = 0 , 

(44) 


where the notation for the above conplings is shorthand for 

The one-loop threshold corrections to Xk{Ms) in the MSSM from box and triangle dia¬ 
grams are tabnlated in, e.g. Ref. [5U] : 




K 


92+9y^ 


12 7 

3hlil^ - 3hlAl - hlAl 




Ai;>A2 = ^ 




^ , 4 -4 ^ 7 4 -4 

2 g 


K- 


.92+9y " 


3h?^ - - Kfl 


2 .A 2 


,2 r.2 


4 

(3 — Al) A 3hl{3 — Al) + h}(3 — Al) 


+ ^hlhl 


3{At + AhY — {(j? — AtAiiY — 6/i^ 


92 - 9y" 


— k^) + ~ k^) + ~ k^) 

Z , ^ L 

^th-^4 = gh^ (3 - Al) + 3hl{3 - Al) + h}(3 - Al) 


K 


- '^hlh^ 


t 


i^9l " 


2 
K 

2 2 L 

K 
6 


3{A, + A,f-{fi^-A,A,f-Qfi^ 


ShiiAf - ii‘) + ihliSl - i,‘) + IliiAi - i^‘) 


3h^A^ + 3/7^74| + 


(45) 

(46) 

(47) 

(48) 


(49) 

(50) 
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At + 3hlAb(Al — 6 ) + h*Ar{A^ — 6 ) 

D L 

^th^'^7 = qA 3hfAt{^ — 6) + 3h^fi^Af, + h^fi^A^ 


(51) 

(52) 


where At = At/Ms, fi = fJ^/Ms, the Yukawas ht^^r are given in Eqs. (25-33), and all 
parameters are in the MS schemej^ 

In addition, there are self-energy corrections to the Higgs bosons which, after redehnition 
of the Higgs fields, give rise to one-loop corrections to the quartic couplings: 

3h‘lfi^ + 3hlAl + 


A ( 1 ) ^ _ f^9l + 9l 
.{1). f^ 9^ + 9Y 

$ 2 g 2 

a(i)x _ ^9I-9y 


3hiAi + 3hifi^ + 

^ht{A^ + jj^) + 3hl{Al + jj^) + h^(^ + /i^) 


3hl{^t + A^) + ^hl{Al + p^) + + /i^) 


Ml). _f^92^ 

eYi 

aL^Ar = aL^Ar = aL^At = 0 


(53) 

(54) 

(55) 

(56) 

(57) 


-^5 — ^<1. ^6 — ^<l> ^7 

We extend these corrections with additional two-loop hfg"^ terms, which can be extracted 
from the corrections to A in the ~ Ms case 


Mhfgi) ^ _ 

^th 


A = IGK^hts^pgU - 2Xt + - -X^ [> , 


7—3 7-4 
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(58) 


and these are matched to the quartic couplings in Eq. (17) by picking out the terms pro¬ 


portional to (c^, sfj, c^s^t, c^pSp, s^^C 0 ) for (Ai, A 2 , A 345 = A 3 + A 4 + A 5 , Ae, A 7 ), respectively: 


2 Ht + , 

3 * 12 * ■' ’ 


A 


- 1 + -4i) , 

- -4? + 

Mtl) 


(59) 

(60) 

(61) 

(62) 

(63) 


Note that there is an asymmetry Aj:j^* ^ \ when Xt —)■ —Xt, however, this is subdominant to 
the asymmetric contribution from the ht threshold in Eq. (28), which leads to log-enhanced 
corrections to the quartic couplings at the two-loop level. 


^ We have not included the small threshold corrections to the quartic couplings from electroweakinos, which 
involve only gY,g 2 , A. They are listed in Ref. [3S], and we estimate they lower mt by about 0.5 GeV. 
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IV. RG EVOLUTION AND DIAGONALIZING THE HIGGS MASS MATRIX 


From the low energy effective theory point of view, the tree-level values of the quartic 
couplings, as well as the threshold corrections enumerated above, should be defined with the 
gauge and Yukawa couplings at the stop mass scale. More precisely, these values define the 
boundary conditions for the RG evolution of the quartic couplings as well as the Yukawa 
couplings from the scale Ms down to the scale of the GP-odd Higgs boson tjia- 

It is clear that if the scale niA is of the order of the weak scale, the GP-even mass matrix 
elements, Eq. 0 . may be computed by evolving all quartic couplings down to the weak 
scale. The GP-even Higgs masses may then be calculated by diagonalizing the mass matrix 
at the weak scale and adding the proper corrections converting the running masses at the 
weak scale to the pole masses. 

On the other hand, if the GP-odd Higgs boson mass is much larger than the weak 
scale, decoupling of the heavy Higgs bosons should be achieved, and the mixing between 
the nonstandard and standard GP-even Higgs boson, Ca-is, should go to zero as l/m\. 
Therefore, for a heavy supersymmetric spectrum, the effective theory is just the SM below 
the scale uia- The lightest GP-even Higgs mass computed in this case should reduce to the 
one previously computed in Ref. [26] . 

Both limits may be appropriately recovered by evolving the quartic couplings up to the 
scale ruA and computing the matrix elements at that scale in the Higgs basis. At scales 
below niA we simply evolve the (1,1) matrix element by the full two-loop SM RG evolution 
of the quartic coupling up to the weak scale. On the other hand, the corrections of the 
off-diagonal elements coming from the evolution from rriA to the weak scale become relevant 
only for large values of niA, for which Ca-i 3 becomes small and therefore irrelevant from the 
phenomenological point of view. Gonsequently, we consider the evolution of this matrix 
element from niA to the weak scale by resumming the dominant top-induced corrections at 
the one-loop level. For similar reasons, for large values of ttia, the radiative corrections to 
the (2, 2) matrix element coming from the running between and the weak scale, which 
depend logarithmically on the GP-odd Higgs mass, become small compared with the tree- 
level value, which depend quadratically on this mass. For low values of rriA < 500 GeV and 
t /3 > 4, we have checked that our results for the Higgs boson masses (mixing angles) differ 
by less than 0.5 GeV (1%) from the ones that would be obtained by evolving all quartic 
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Observable 


Value 


SU{3)c MS gauge coupling (5 flavours) as{Mz) = 0.1184 ± 0.0007 


Fermi constant from muon decay 
Top-quark pole mass 
W boson pole mass 
Z boson pole mass 
Higgs pole mass 


V = = 246.21971 ± 0.00006 GeV 

Mt = 173.34 ± 0.76 ± 0.3 GeV 
Mw = 80.384 ± 0.014 GeV 
Mz = 91.1876 ± 0.0021 GeV 
Mh = 125.09 ± 0.21 ± 0.11 GeV 


TABLE I. SM observables, collected in Table 2 of j37| . 


couplings until the weak scale. For larger values of itia and lower values of tp, however, 
the effects of decoupling of the heavy Higgs bosons become relevant and for tua 3> 1 TeV, 
the lightest GP-even Higgs boson mass can become signihcantly different from the one that 
would be obtained without decoupling the heavy Higgs bosons. 

As a starting point for the evolution of the quartic couplings, one specihes the SM values 
of the gauge couplings and Yukawa couplings at some scale. We work in the third-generation 
approximation, so the six couplings 93 , 92 , gi,yt,yb,yT are relevant. We use the low energy 
parameters 9 i,yj at the scale of the top-quark pole mass AA, which are extracted from the 
SM observables in Table|^ and have values given in Table [lT|P1 These couplings, along with an 
initial valu^ |37l [391 - H8] of A 0.25, are evolved to the intermediate scale rriA using three- 
loop SM RG equations for 93 , 92 ,9i,yt, X and two-loop SM RG equations for yb,yT- There 
are additional loop contributions to 91 , 92 , yt, X from electroweakinos if p. Mi, M 2 < rriA- [IS] 
Due to their weak couplings, we have only included the dominant one-loop log-enhanced 
contributions from RG running using tree-level gauge couplings of the electroweakinos to 
the Higgs bosons. 

Above the scale rriA, the effective theory is the THDM, and two-loop type H THDM RG 
equations are employed in the running between tua and Ms- These are listed in Appendix |A[ 
and can also be found in Ref. [5^. As above, we have included one-loop contributions to 
the running of 9 i, 92 ,ht,Xk from electroweakinos if < p, Mi,M 2 < Ms- We note that 
for perturbative consistency of the RG running, three-loop RG equations should be used; 
^ Unlike in Refs. [5^[2H], we use the NNLO value of yt{Mt) instead of the NNLO-|-N^LO QCD value because 

we use three-loop SM RG equations below but only two-loop THDM RG equations above ruA- 
^ We have checked that the final values for Mh do not have a strong dependence on the initial condition 

for \{Mt) if it is chosen to correspond to a value of mh{Mt) = Au^ ~ 100-150 GeV. 
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however, these are not known for the THDM. Also, inclusion of the three-loop order RG 
equations in the SM running has a small effect, and we expect the same holds for the THDM. 
To determine the approximate values of the MSSM gauge and Yukawa couplings at the high 
scale Ms, we run the couplings up to Ms setting = 0 (A; = 1,..., 7) in the running; this 
has subpercent level effects on the running of the Yukawas and the gauge couplings. At Ms, 


we calculate the threshold corrections to the Yukawas according to Eqs. (25)-(33), and use 
these results in the expressions for the MSSM values of the in Eqs. (41)-(57, (59)-(63)|^ 
With these values of Xk{Ms), we use the full type II THDM RG equations in the running 
back down to m^. The matrix elements of Ai'jj in Eqs. ( |l^ [T7]-[T^ are computed, and the 
value of gn is used as the boundary value for X{mA) for the SM RG running from tua to 
Mt- Ml'jj is then diagonalized at Mt, and the running masses and the mixing angle 

(5 — a are computed. 

The running mass rrih is converted to the pole mass M^ using the SM one-loop formula 
as in Ref. [20], in which SM MS running couplings are used. 


Ml = X{Mt)v\Mt) + 3yl{4:mf - ml)Bo{mt,mt,mh) - ^Am^ 


2 - ^ 


1 


V 

~4 


2 r 


3g2 - 4Xgl -F 4A^ BQ{mw, mw, mu) + -glv'^ 


9l 


A log 


m 


w 




.2 r 


(64) 


3(92+Sr) - 4A(92 + Sr) + 4A Ba{mz,mz,mk) 


nir 


+ 4 (^2 + ai-y [{9i + - A (log ^ - 1 

where Bq is the one-loop Passarino-Veltman integral 


Q^=M^ 


Bo{mi,m2,m3) = 


log 


2 -, 


(1 — x)mi + xml ~ ^(1 ~ x)ml 

0^ 


(65) 


We have checked that this gives similar results as the two-loop conversion using parameter 
values in the on-shell (OS) scheme, as in Ref. [37]. We have not included contributions from 
light electroweakinos in the conversion formula, but we expect these contributions to be 
subdominant to those from the SM. For low values of mA and and large values of Ms, the 
top Yukawa yt can deviate sizably from the coupling of the physical Higgs to the top, ghtt in 


Eq. (35); however, we checked that the shift in M^ when substituting ghu for yt in Eq. (64) 


We use tp{mA) as an input, and run vi and V 2 to Ms using two-loop RG equations for the anomalous 
dimensions. When calculating the MSSM values at Ms, we use tp{Ms) = V 2 {Ms)/vx{Ms)■ 
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9i 

gi{Mt) 

Vj 

yj{Mt) 

93 

1.1666 

yt 

0.94018 

92 

0.64779 

Vb 

0.0156 

gy = y/ 3755 'i 

0.35830 

Vt 

0.0100 


TABLE II. Values of SM parameters at Q = Mt computed in the MS scheme. g 2 ,gY, and yt are 
computed at NNLO, and the SU{5) normalization relates <71 to the SM hypercharge conpling gy- 
The value of g^ is obtained using three-loop QCD matching to the SM. We have nsed the two-loop, 
five-flavour MS RG equations in the broken phase from j38| to rnn rur from their initial values 
mb{mb) = 4.18 GeV, M,- = 1.777 GeV [51]. For more details, see [37] . 

is less than 0.5 GeV. Similarly, in the scenarios we investigated, we expect the difference 
between the rnnning mass and the pole mass of the heavy Higgs to be small dne to its larger 


mass and its rednced conpling of the heavy Higgs to the top qnark in Eq. (36) 


It is instrnctive to consider the dominant one-loop contribntions to the GP-even Higgs 
matrix elements in the Higgs basis. This was discnssed in detail in Ref. [3l] (for the CP- 
violating case see Ref. [36]), in which it was shown that 


Sv^sih'f 

giiv = m\clp + 


Svr^ 


gi2V^ = -S 2/3 rn\c2^ 


In ( ^ ) + 

mi 


1 -^ 
12 




In 


, XtiXt + Yt) 


\mi 


+ 


X^tYt 


12 


( 66 ) 


(67) 


where Xt = Xt/Ms, Xf = At — is the stop mixing parameter associated with the 
conpling of the SM-like Higgs to the stops, 1) = Yt/Ms and Yt = At + ntjs. 

From the above, the mixing angle C/ 3 _q may be compnted, 

_ -gi2v‘^ 

^0—a — 


- ml){mij - guv'^) 
For valnes of mu larger than the weak scale, one can show that 


( 68 ) 


— 


-1 


my-ml 


2 2 ^?>mtXt{Yt-Xt) 

rrih - mzC2p + - 



(69) 


47^2^2 

and therefore, all dominant radiative corrections to the mixing angle, which come from the 
renormalization of the qnartic conpling A 2 at scales above mA, may be absorbed into the 
dehnition of the Higgs mass mh- The remaining terms are proportional to /tXttan/3, vanish 
for maximal mixing Xt = 6, and cannot be absorbed into a redehnition of mh- 
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A. The hMSSM scenario 


The fact that the dominant corrections to the lightest CP-even Higgs mass and the 
mixing parameter have a common origin motivated the authors of Ref. [S2] to dehne the 
“hMSSM scenario”, in which one assumes that for a given value of and ruA the proper 
Higgs boson mass may be obtained by choosing appropriate stop masses and mixings. As 
discussed above, the dominant radiative corrections to the Higgs boson mass and the CP- 
even Higgs mixing may be absorbed into the definition of the Higgs boson mass irth. More 
precisely, in the hMSSM scenario, only radiative corrections to the quartic coupling A 2 are 
considered, namely Ai = —A 345 = A 2 = -|- AAd 22 /('^^' 5 ^) Ag = Ay = 0 . 

One can easily show that, in order to achieve the proper Higgs pole mass Mh, the radiative 
corrections must be given by 




Ml(ml + M| - Ml) - mlMlclg 

Af|4 + - Ml 


(70) 


The heavy CP-even Higgs mass is given by 

2 _ + m\sl) - 

- Mlcl + m\sl - Ml 


(71) 


Once these expressions are considered, the CP-even Higgs mixing angle is given by 


a = — arctan 


(M| + m\)cfiSp 


Mlcl + m\s^ 






Ml 


(72) 


It is straightforward to show that, for values of ijia larger than the weak scale, the mixing 


angle in this approximation agrees with the one presented in Eq. (69), in which the last term 
inside the bracket is neglected. 


From Eq. (69), one can then identify the main difference between our approach and the 
hMSSM approximation. For low values of /i, the main difference is associated with the 
radiative corrections to the quartic couplings Ai and A 345 . For moderate values of tp, the 
main logarithmic corrections to these couplings are governed by weak couplings, and hence 
are small compared to the dominant corrections absorbed into m/j. 

On the other hand, for sizable values of fi and moderate values of tp, the last term 


in Eq. (69) may become relevant and therefore we expect the hMSSM scenario to fail to 
accurately describe the Higgs phenomenology in this case. The difference will be maximal 
for sizable values of fi and Xt away from the maximal mixing value. Beyond the difference 


16 








associated with the /i-induced radiative corrections, the liMSSM works under the assumption 
that one can always choose the supersymmetry-breaking parameters so that the proper Higgs 
mass is obtained. As we shall show, this is not always possible: for low values of ttia and 
t/ 3 , there are regions of parameter space for which obtaining the right Higgs mass would 
demand pushing the supersymmetric particle masses close or above the Planck scale, where 
the effective low energy supersymmetry description is no longer valid. 

We can also examine the Hhh coupling. We dehne the Feynman rule for the vertex as 
iguhh- This is given by |3ll |53] 


dHhh — 


-Sns^c^l — -\- SapCajS (Ai — A345)sQ,y3 -|- tji{^Q{2Cal3 + Sap) — X^{2Saj3 + Cq,^) 

+ — + A345) -f — X'iibVCa-p , 

(73) 


where Sap = (—Sa/c/ 3 ) and Cajs = both of which tend to 1 in the alignment limit. 

The above expressions can be compared to the expression given in the hMSSM approxi¬ 
mation 


gUhh — 


Ml 


2S2q:S/ 3+Q! C2Q:Cy3^Q -|- 3 


AAf 


22 2 


Ml sp 


(74) 


which can be recovered from Eq. (73) when the radiative corrections to A 2 alone are con¬ 
sidered. Hence, as with the mixing angle a, we expect the hMSSM to provide a better 
approximation to the correct results provided fi is small. 


V. NUMERICAL RESULTS 

The results of our analysis are presented in Figs. iii and 1^ In Fig. we present 
contour plots of the lightest CP-even Higgs mass for rriA = 200 GeV in the Ms-tp plane 
for small values of the stop mixing parameter, Xt = 0 , and for values close to maximal 
mixing, Xt = a/O- In addition, we compare values of obtained for jj, = Ms-, for which 
the chargino and neutralino contributions to the Higgs mass decouple below the scale M 5 , 
with the ones for low values of /i = 200 GeV, for which the corrections to the Higgs mass 
induced by RG-evolution effects of charginos and neutralinos become relevant. We see that 
in order to obtain the proper value of the Higgs mass at low values of tp ~ 2, low values of 
fi of the order of the weak scale and large values of Ms of the order of the GUT scale are 
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Mi„mA=200 GeV,ii=Ms,X,=0 


Ml,, m^=200 GeV, ti=200 GeV, X,=0 
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Log,„[Ms/GeV| 
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FIG. 2. Contour plots for in the plane tjs, Ms with mA = 200 GeV, Mi = M 2 = ^, Ai, = = 

At- The outer black, dot-dashed lines are contours of Mh = 122,128 GeV as labelled. The blue, 
dashed lines correspond to = 124,126 GeV, and the central black, solid line to = 125 GeV. 
Plots in the top (bottom) row have Vt = 0 (\/6), and plots in the left (right) column have /x = Ms 
(200 GeV). 


necessary. We also note that for tp < 1.5, values of Mh = 122 GeV may not be obtained 
even if the supersymmetric spectrum is pushed to the GUT scale. 

The values of the Higgs mass at tua = 200 GeV are heavily susceptible to Higgs mixing 
effects. In contrast, we show in Fig. [^contour plots of the lightest GP-even Higgs mass for 
tha = 300 GeV and similar supersymmetry breaking parameters as in Fig. The qualitative 
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FIG. 3. As in Fig. with niA = 300 GeV. 


behavior is the same as in the previous case, but the proper Higgs mass is achieved at lower 
values of Ms- In particular, for low values of p, values of = 1 no longer demand sparticles 
above the GUT scale, a result that is independent of the stop mixing parameter. 

It is relevant to show the previous results for values of Ms of order of the TeV scale, as 
expected if supersymmetry is related to the mechanism of electroweak symmetry breaking. 
The results are presented in Figs. |^and[^ in which the values of Ms are restricted to vary 
between 1 TeV and 30 TeV. For tua = 200 GeV, it is clear that Mh = 125 GeV cannot 
be achieved with values of tp < 4; similarly, requiring values of Ms of the order of 1 TeV 
demands large tp and values of Xt close to the maximal mixing values. As expected, the 
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FIG. 4. As in Fig. but with Ms restricted to the 1-30 TeV range and modified ranges in 
Shading has been added between the contours for visual clarity. 


values oitp necessary to achieve the proper Higgs mass increase for lower values of and 
large values of /r, due to mixing and chargino and neutralino effects, respectively. 

In Fig. 1^ we plot as a function of to show the effect of mixing in the Higgs mass 
matrix at different values of rriA, t^, hxing p = 200 GeV. The different curves correspond to 
choice of Ms between 1 and 10 TeV. At low = 5, the effect of mixing for ttia = 200 GeV 
is pronounced; the value of Mh with ttia = 200 GeV is between 2-3 GeV lower than with 
niA = 500 GeV, which approximates the decoupling limit. For higher values of = 20, the 
difference between the respective curves for the two values of ruA falls to less than 0.5 GeV. 
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Ms [GeV| Ms [GeV| 


FIG. 5. As in Fig. but with Ms restricted to the 1-30 TeV range and modified ranges in 
Shading has been added between the contours for visual clarity. 

The logarithmic dependence of on Ms is evident in these plots: increasing Ms from 1 to 
5 TeV increases M^ by approximately 12 GeV, while doubling Ms from 5 to 10 TeV yields 
a more modest 3-4 GeV change. We also note that the maximum M^ achieved for M 5 = 1 
TeV is Mh ~ 126.1 GeV in the lower right panel. Within uncertainties, this agrees with 
results previously found in the niA = Ms case in Ref. [26] . 

In Fig. we have plotted Mh in the high-scale SUSY scenario, with large = 20. For 
Ms = 2 TeV, the dashed blue curve, we obtain Mh = 126.5 GeV. Also, in contrast to 
the bottom-right panel of Fig. Mh = 125 GeV is no longer achieved for Ms = 1 TeV 
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m^=200 GeV, tan/3=5, ^=200 GeV 



mA=200 GeV, tan^=20, /^=200 GeV 



Xt 

-Ms=1 TeV . Ms=2 TeV 

. Afs=5TeV - Ms=10TeV 


mA=500 GeV, tan^=5, fi=200 GeV 



mA=500 GeV, tan^=20, /i=200 GeV 



X, 

- Ms=^ TeV . Ms=2 TeV 

. Ms=5TeV - /lfs=10TeV 


FIG. 6. Mfi vs Xt for uia = (200,500) GeV in the (left, right) columns, tjs = (2,20) in the (top, 
bottom) rows, Ai, = Ar = Ms, and fj, = Mi = M2 = 200 GeV. The four curves are for Ms values of 
1, 2, 5,10 TeV from bottom to top. The vertical grey dashed line indicates the value at the one-loop 
maximal mixing value Xt = \/6- The horizontal light grey box is the la band M^ = 125.09 ± 0.24 
GeV. 

at maximal mixing without light electroweakinos. We can compare with the recent results 
produced by the SusyHD code of Ref. [2H]. Our values are < 1 GeV higher than the central 
result of Ref. [2H]. Part of this discrepancy is attributed to the use of the lower value of 
yt{Mt): if we instead use the NNLO + N^LO QCD value yt^N^LO QCD(^t) = 0.93690, Mh is 
lowered by 0.5 GeV. The remaining small difference may be explained by the more complete 
calculation of thresholds in the ~ Ms case of Refs. |2ni EH]- 


VI. COMPARISON TO PREVIOUS RESULTS 

In this section, we compare our results with the results obtained in the liMSSM scenario 
as well in the FeynHiggs version 2.10.2, in which relevant logarithmic effects to the SM 
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mA=Ms, taiv5=20, fi=Ms 



X, 


- Ms=^ TeV . Ms=2 TeV 

. /Ws=5TeV . /lfs=10TeV 


FIG. 7. As in Fig. with ruA = Ms, tp = 20, = Ms, and ^ = Mi = M 2 = Ms- 

quartic couplings are resummed in order to increase the accuracy of the results at large 
values of Ms [MIES]. 

In Fig. we present the comparison of our results with the hMSSM approximation for 
sizable values of /t = 2 and values of Xf = —1.5 and Xf = 2.8, away from maximal mixing, 
for which the hMSSM results are expected to show a worse approximation to the correct 
results than for low values of p at moderate or large values of The results of our compu¬ 
tation for the mixing angle a and the heavy GP-even Higgs mass are presented in the left 
and right panels with red dotted lines, while the blue lines represent the relative and abso¬ 
lute differences of these quantities with the ones computed in the hMSSM approximation. 
We present our results for Ms = 5 TeV, for which the correct values of the Higgs mass, 
represented by black solid, dashed and dotted lines, may only be obtained for moderate to 
large values of tp in this region of parameters. Differences in a of the order of 10%-20% 
are obtained for moderate values of tp and values of the heavy GP-even Higgs bosons of 
the order of the weak scale. Since the mixing angle controls the coupling of the lightest 
GP-even Higgs boson to fermions and gauge bosons, relevant modihcations of the Higgs 
phenomenology are expected in this region of parameters. Similarly, the heavy GP-even 
Higgs boson mass may be affected by values of a few to 10 GeV in this region of parameters. 

In Fig. we present in the upper panels similar results but for Xt = 2.8 and large values 
of Ms = 100 TeV for which lower values of ~ 4 are required to obtain the correct Higgs 
masses. We see that in this case, in the relevant region of parameters, the agreement is 
improved compared to the large tp case, with differences in a of the order of a few percent 
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niA [GeV] aiA [GcV| 




FIG. 8. Difference between hMSSM and effective THDM calculations of a and Mh in the plane 
tjs,mA for Ms = 5 TeV, Xt = [2.8, —1.5] [top, bottom], /r = Mi = M 2 = 2 M 5 , and Ai, = A^- = At. 
From the bottom to the top of each plot, the light grey lines (dot-dashed, dashed, and solid) 
correspond to Mh = (122,124,125,126,128) GeV. Red, dashed lines in the plots in the left (right) 
column are contours of a [Mh) computed using the effective THDM. Solid and dashed blue lines 
in the plots in the left (right) column are contours of (a —ahMSSM)/|«| {Mh — Mu^^mssMi in GeV); 
dashed lines indicate values halfway between adjacent solid lines. 


and differences in m/f of the order of a few GeV. In the lower panels, we present results for 
lower values of fi and Ms = 10 TeV, for which values of t /3 ~ 5 lead to the proper Higgs 
boson masses. We see that due to the smaller values of (1 and tp, the differences with the 
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Aff/|Q-|,Mj=10TeV,X,=2,//=l TeV 




„U [GeVI 

Am//,Mj=10TeV,A',=2,^=1 TeV 



FIG. 9. As in Fig. with Ms = 100 TeV, Xt = 2.8, jl = 2 [top] and Ms = 10 TeV, Xt = 2, ^ = 1 
TeV [bottom]. 


hMSSM reduce to values of at most l%-2% in this case. 

A “low-tan/3-high” scenario, in the region of 1 < < 10, 150 GeV < tua ^ 500 GeV, 

has been presented by the LHG Gross Section Working Group [M] with values for a subset 
of the MSSM parameters necessary to achieve = 122-128 GeV in FeynHiggs. In this 
scenario, a simplihed, heavy MSSM spectrum above the scale Ms is assumed; other MSSM 
parameters are chosen as Af = 2 TeV (/ = b,T,c, s, ^,u,d,e)] M 3 = Ms; M 2 = 2 TeV; 
/X = 1.5 TeV; and Mi = M 2 ■ |tan^ 6 *w ~ 950 GeV fixed by the GUT relation. M°® and 
the values of the stop masses and mixing parameters in the on-shell scheme used in 
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FeynHiggs, are then chosen to achieve Mh in the desired range. 

We have used one-loop conversion formulas [HIES] to change in the OS scheme 

to Ms,Xt in the MS scheme, which are the parameters used in our calculation. The maxi¬ 
mum value specihed is 100 TeV, which is used for points with low tua ^ 200-250 GeV 
and tjs < 1-3. Maximal mixing is chosen for points in the region tp < 2. In FeynHiggs, this 
corresponds to the choice = 2Mg^] in the MS scheme, the output value of Xt should 
be close to the maximal mixing value for which as a function of Xt achieves its 


maximum (e.g., in Fig. 


X 


lies close to the one-loop value = \/6Ms). For 


t,max 


selection of these points, we performed a similar scan and found that the output values of 
Xt yield Mh{Xt) values within 0.5 GeV of the maximal mixing values. 


Our results in this scenario are shown in Fig. IT The top-left panel shows that across the 
range of parameter space using the tabulated values of Ms and Xt, Mh ^ 123 GeV using the 
effective THDM calculation. In the top-right panel, the discrepancy between our calculated 
value of Mh and that of FeynHiggs is clearly exhibited: for much of the parameter space 
above tp ~ 6.5, our calculation of Mh is about 2 GeV lower. Between 5 {tp ~ 2—4), 

this disagreement worsens to 3-5 (5-10) GeV. This can also be seen in the lower-right panel 
of Fig. 1^ where for niA = 200 GeV and lower values of /x = Mi = M 2 = 200 GeV, Mh = 122 
GeV is not achieved for Mg = 100 TeV until tg ~ 3. The effective THDM calculation yields 
a higher value of the Higgs mixing angle a compared with FeynHiggs, but the two are 
in agreement at the level of 5% except for a region < 300 GeV and < 6. Below 
~ 3 and tua ~ 225 GeV, the fractional difference reaches 10%-12%. The values of the 
heavy Higgs mass tuh are only signihcantly discrepant, more than 5 GeV, for low < 2.5, 
although for tp < 1.5, rriA ^ 250 GeV, our calculation of mn is more than 10 GeV lower. 


We can estimate how much of the differences in Fig. [^are due to the use of a different 
boundary value for the top Yukawa yt{Mt), for which FeynHiggs uses the one-loop SM 
MS frit running valuej^ In Fig. 


11 


we reproduce the results in Fig. 10, except that we use 


yt,T<!ho{Mt) = 0.95113 as the boundary value for the RG running. In the top-right panel, Mh 
in the region above tp ~ 5.5 (4.5) now agrees to within 1 (2) GeV; however, discrepancies 
larger than 5 (10) GeV still exist for tp 3 (2). Likewise, there are modest reductions in 
the differences in a and rriH across the parameter space. The remaining differences between 
® For consistency at two-loop order, only the one-loop terms involving gj,,yt are employed in FeynHiggs 
to obtain NLo(Mt) = 0.962. 
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FIG. 10. Top row: Density plots for Mh calculated using the effective THDM (left) and the 
difference between the left plot and the calculation of Mh using FeynHiggs, for the low-tan/3- 
high scenario (right). From top to bottom, the (dotted, dashed, solid, dot-dashed) black curves 
correspond to differences of —(2,3,5,10) GeV, respectively. Bottom row: The difference in mn 
(left) and fractional difference in a (right) calculated using the effective THDM and FeynHiggs. 
In the left plot, from the upper right to the lower left, the (dashed, dotted, dotted, dashed, solid, 
dot-dashed, dot-dashed) black curves correspond to differences of (2,1, —1, —2, —5, —10, —15) GeV, 
respectively. In the right plot, from top to bottom, the (dotted, dashed, solid, dot-dashed) black 
curves correspond to differences of (1,2,5,10)%. 


the FeynHiggs results and our results could be explained by the different resummation 
method implemented in FeynHiggs in which the THDM effects are ignored. 


We turn now to the comparison with the hMSSM in this scenario, shown in Fig. 1^ We 
use Eqs. (|7g-(|7|), inserting the value of Mh obtained from the effective THDM calculation. 
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FIG. 11. As in Fig. 10 except using the boundary value yt,NLo(Alt) = 0.95113 for the RG evolution. 
In the top right plot, from top to bottom, the (dotted, dotted, dashed, solid, dot-dashed) black 
curves correspond here to differences of —(1, 2, 3, 5,10) GeV, respectively. 


The fractional difference in a between our calculation and the hMSSM is less than 4% 
between the two calculations. Likewise, there is minimal deviation in m^, except in the 
small corner of parameter space at f/? ~ 1, < 200 GeV, where the disagreement reaches 


the 5% level. As was discussed in Sec. IV A sizable values of /r are needed for the hMSSM 
approximation to break down; however, throughout the parameter space of the low-tan/3- 
high scenario, Mg- Finally, we note that if instead the value of Mh from FeynHiggs 
is used in the hMSSM equations, we see a similar level of disagreement between the hMSSM 


and our calculation as in Fig. 10 


We can also test the formulas for the gnhh coupling, Eqs. (73)-(74), in the low-tan /3-high 


scenario. In Fig. 13, we show the results of our calculation and the fractional difference 
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FIG. 12. Density plots for the difference in mn (left) and the fractional difference in a (right) cal¬ 
culated using the effective THDM and the hMSSM approximation, for the low-tan /3-high scenario. 
In the left plot, from the upper right to the lower left, the (dashed, dotted, dotted, dashed, solid, 
dot-dashed) black curves correspond to differences of (2,1, —1, —2, —5, —10) GeV, respectively. In 
the right plot, from the upper right to the lower left, the (dotted, dashed, solid, dot-dashed) black 
curves correspond to differences of (0,1,2,3)%, respectively. 

with the hMSSM using the effective THDM value of Mh- Fractional deviations of less than 
6%-7% are observed. As above, differences between our calculation and the hMSSM when 
the FeynHiggs value of Mh is used reach 30% at low tan/3 and larger values of tua- 

The dominant SM uncertainties come from the inputs yt, at Mj. The uncertainty from 
as{Mt) is subdominant as it enters at two-loop order for Mh in both the RG running of 
yt, Xi and in the threshold contributions. The uncertainty from yt has two sources: one 
from the experimental measurement of the top-quark pole mass Mt, and the other from the 
conversion of Mt to the MS top Yukawa yt{Mt). An estimate of the uncertainty from the 
value of Mt can be found in the ~ Ms case in [2^, where it was shown that using the la 
high and low values of Mt shift Mh by about 1 GeV. As previously discussed, the use of the 
NLO, NNLO, or NNLO+N^LO QGD values of yt{Mt) can shift Mh by 1-2 GeV. There are 
also uncertainties from varying the renormalization scale in the effective potential, from 
subleading two-loop threshold corrections to A^, and from higher-dimensional operators, but 
we expect these contributions are subdominant to those from the SM. For a more detailed 
discussion of uncertainties, see Ref. [28] . 
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FIG. 13. Contour plots of Xnhh = QHhh/v computed in the THDM [left] and the fractional 
difference between Xnhh computed in the THDM and the hMSSM [right]. 

VII. CONCLUSION 

In this article, we have computed the mass and couplings of the lightest UP-even Higgs 
boson in the MSSM, considering large values of the masses of the scalar quarks, and inter¬ 
mediate values of the UP-odd Higgs mass. We performed these calculations using effective 
theory techniques and resumming the large logs appearing above and below the UP-odd 
Higgs mass scale. We worked in the Higgs basis and showed that provided the threshold 
corrections to the off-diagonal UP-even Higgs mass matrix element are small, all relevant 
radiative corrections may be effectively absorbed into the definition of the lightest UP-even 
Higgs mass. This situation occurs for moderate or small values of the Higgsino mass pa¬ 
rameter fi and/or of the trilinear stop mass parameter At, and the resulting UP-even Higgs 
boson masses are well approximated by the hMSSM scenario. On the other hand, for sizable 
values of /i and At, the alignment condition may be realized, in which case our results differ 
signihcantly from those in the hMSSM method. 

The Higgs masses computed in our work tend to be lower than the results obtained by 
FeynHiggs, which implements a different resummation method, and may differ by a few 
GeV or more. The difference may be traced to our use of an effective THDM theory at 
scales above rriA and also a higher-order computation of the relation between the running 
and the on-shell top-quark mass. 


30 




















Our calculation of Mh leads to lower bounds on tp for low values of tua in order to achieve 
Mh = 125 GeV: for = 200 GeV, we find > 3.4 (2.0) for /i = Ms (200 GeV). These 
bounds are due to the appearance of large mixing effects that push the lightest GP-even 
Higgs mass down and cannot be overcome by the positive contributions to Mh from both the 
stop threshold corrections (even with stops as heavy as Mqut) and from radiative corrections 
from light charginos and neutralinos. We also note that for values of f/? ~ 1, a Landau pole 
of the top-quark Yukawa coupling may be induced at low values of niA- 

Finally, we would like to stress that our work has been restricted to the computation 
of the neutral Higgs masses in the MSSM, without taking into account any experimental 
constraints beyond the measured value of the Higgs mass. While constraints from precision 
measurements, new physics searches, Higgs, dark matter, and flavour physics will lead to 
relevant bounds on the values of the free parameters of the theory, in this article we have 
focused on the Higgs mass computation for arbitrary values of those parameters. Our results 
should be complemented with a careful analysis of the experimental constraints and can also 
be used to determine in a more precise way the bounds on the free parameters of the model 
coming from those constraints. We reserve this analysis for a future publication. 
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Appendix A: Two-loop RGB’s in the Type II THDM 


The /9-functions for a coupling are 


, OO 


(Al) 


n=l 


where t = logQ with Q the renormalization scale, k = l/(47r)^ is the loop factor, and 


is the nth loop ,9-function for g. We have extracted these equations from the program 


SARAH, version 4.2. |57j Below, we list the two-loop RG equations for the type II THDM 
in the third generation approximation. Ng is the number of fermion generations and 6x 
is the Heaviside function for the mass X. These equations were also listed in Ref. [SU] . 
with which we hnd minor differences; we use different conventions for three parameters 
Al = 2Ai, A 2 = 2 X 2 , gf = 5g'‘^/3, where g', Ai, A 2 appear in Ref. [50] . 

1. Gauge couplings 

Hypercharge coupling gi in the SU{5) normalization, with gy = ^gf'- 



(A2) 


(A3) 


Weak gauge coupling g 2 '. 



(A4) 


(A5) 


Strong gauge coupling g^: 



(A6) 


(A7) 
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2. Yukawa couplings 


Top Yukawa ht'. 
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Bottom Yukawa hb'- 
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Tau Yukawa hr- 
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3. Anomalous dimensions 


Down-type Higgs 71 = d\ogVi/dt\ 
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Up-type Higgs 71 = d\ogV 2 /dt: 
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4. Quartic couplings 


= l 2c/^ + (^gl + - 12ht - Ah^ 

+ 12 A 2 + 2(A3 + A 4 )" + 2 A 2 + 2\l + 2A\l - 4Ai7i 
— ^glOM 2 + (^fji — 4 Ai 7 {^^ 

( 2 ) 291 6 303 4 2 1719 2 4 3537 ^ 

/^Ai g 92 4Q 929 i 200 1000^^ 

- —5^2 12/i^ + 4 /i^ + 17 Ai — 2 o( 2 \^ + A4') + ~ 5 ' 25 'i 36 /i^ + 44 /i^ + 39 Ai + 2OA4 

o L \ / J 2(j L 

+ 200 ^^ 20/i^ — lOO/i^ + 217Ai + 20^2A3 + A4^ 

— IQg^h^ Ah'^ — 5Ai + — 5^2 5Ai ^3/i^ + + 4^9A^ + 4 A 3 + 4 A 3 A 4 + A 4 + 18Ag^ 

+ + 25Ai^ + 3h‘^ ^ — 16/i^ + 25Ai^ 

+ 12 ^ 9 A 4 + 4A3 + 4A3A4 + 2A4 — A5 + 18 Ag^ 

- 12/1^ [2A^ + 2A3A4 + A^ + A^ + 6A^j + Sh^^hl Uhl - 3 Ai 

-3hl - 20ht + hlXi + 2AXl + 2AXl -hi - 20hl + hlXi + 2AXl + 2AXl 

- 2 X 1 [ 39 A? + lOA^ + IOA 3 A 4 + 6X1 + 7 A 5 + 159A^ - 3X1 

- 4 A 3 Ua^ + 6 A 3 A 4 + 8 A^ + lOA^ + 33A^ + 18 AgA 7 + 9A? 

- 4 A 4 [3A^ + llA^ + 35A^ + I 4 A 6 A 7 + 7A^j - 4 A 5 [37A^ + 10AgA7 + 5A? 

(A 19 ) 


35 



I^\2 “I ‘^92+ {^2 + “ 12 /lj 

+ 12A2 + 2A2 + 2(A3 + A 4 )" + 2A2 + 24A2 - 4 A 272 
6 9 1 

~ %2^M2 + gfi'2fl'l^M2^Mi + ~ 4A272^^ 

( 2 ) _ 291 6 303 4 2 1719 2 4 3537 ^ 

/^A 2 g 92 4 Q 929i 200 1000^^ 

- ~92 12^i + I 7 A 2 — 20(^2A3 + A 4 'j + ~5'2fi'i 34 / 1 ^ + 39 A 2 + 2 OA 4 

o L \ / J 2(j L 

- [re/i? - 217 A 2 - 2o(2A3 + A 4 ) 

— IQg^h^ 4/ij — 5 A 2 + “(72 15/i^A2 + 4^9A2 + 4 A 3 + 4 A 3 A 4 + A 4 + ISA^^ 

+ ~ hf (32hf — 85 A 2 j + 12 ^9A2 + 4 A 3 + 4 A 3 A 4 + 2 A 4 — Ag + I 8 A 7 j (A 21 ) 

+ ?>hl [ 20 / 1 ^ - hl\2 - 24A^ - 24A^j + ?>hlhl Uh] - 3 A 2 
- ^ 12 /ig + 2 A 3 + 2 A 3 A 4 + A 4 + Ag + 6 A 7 

- 2 A 2 [39A^ + lOA^ + IOA 3 A 4 + 6 A^ + 7A^ - 3A^ + 159A? 

- 4 A 3 Ua^ + 6 A 3 A 4 + 8 A^ + lOA^ + 9A^ + I 8 A 6 A 7 + 33A^ 

- 4 A 4 [3A^ + llA^ + 7A^ + UAgAy + 35A?j - 4 A 5 [sA^ + lOAgAy + 37A^ 


36 



/3i‘’ = ^[29j+(9^-f9?)] -12« 

+ 2(Ai + A2)(3A3 + A4) + 4A3 + 2A4 + 2A5 + 4A0 + I6A6A7 + 4Ay — 2A3 7i ^ 4“ 72 ^ 

- 5c/^0M2 - ^gjglOM^OMi + Of, - 2A3 + 72^^^ 

(A22) 

^(2) 291 0 33 4 2 909 2 4 3537 0 

l3[ = -On 4-5'2fl'? 4-Oofi'i-S'? 

/^As g y 2 40^2^1 200^^^^ 1000^^ 

~ ~^g2 4" 3/i^ + — 15 ^Ai + A 2 ^ 4- “^-^3 ~ IOA 4 

~ YQ^ 2 fi'i 21/it 4- 9/i^ + 11/i^ + 5^Ai + A 2 j —^A 3 + 6 A 4 

+ 2 ggfi'i — 38/it 4“ 10/i^ — 50/i^ + 30 ^Ai + A 2 ^ 4- I 97 A 3 + 2 OA 4 

- 8 c /3 Sh^hl - 5A 3 (^/i? + hfj 

4" %2 g-^3^3/it + Shf + /i^^ + 3^Ai + A2^ ^2A3 + A4^ + ^A3 — A4^ + I8A0A7 
+ gfi'i - ^hfhl + -A3 ^17/it + 5hl + 15/i^j 

+ 12 ^Ai + A2^ ^3A3 + A4^ + 6 ^A^ — A4 + 2A0 + A0 + I6A0A7 + A^^ 

- ^^3 3/it + 3/ife + + /it/ife 36 ^/it + hfj + I5A3 

- 6/it 6A2A3 + 2A2A4 + 2A3 + A4 + A5 + 8A6A7 + 4A7 

- ^6/i^ + 2/i^^ 6A1A3 + 2A1A4 + 2A3 + A4 + A5 + 4A0 + 8A6A7 

- ^At + A 2 ^ I 5 A 3 + 4 A 4 — 2^Ai + A 2 ^ I 8 A 3 + 8 A 3 A 4 + 7 A 4 + 9 A 5 

- 2Ai [31A^ + 22A6A7 + llA^j - 2A2 [llA^ + 22A6A7 + 31A? 

- A3 [l2A^ + 4A3A4 + 16A^ + 18A^ + 60A^ + 176A6A7 + 60A? 

- A4 [l2A^ + 44A^ + 68A^ + 88A6A7 + 68A?j - A5 [68A^ + 72A6A7 + 68A? 

(A23) 


37 



A 4 : 


— ~^ 929 \ 2 A 4 (Ai + A 2 + 4 A 3 + 2 A 4 ) + 8 A 5 + lOAg + 4 A 6 A 7 + lOA^ 

0 


2 A 4 


7!n7?^' 


4 12 2 2 

^92^M2 ^929i^M2^Mx 

5 


- 2 A 4 


7i,x + 72,X 


(A24) 


A? = -^9t9t - ^9l9t + A 4 


42 


657 


231 4 1413 4 


-^9l + 


200 


9i 


3 

To 


51 


929i \ 42 /i^ + 18/^5 + 22 h‘^ + lOAi + IOA2 + 4A3 H— 


+ Sfl'a 


+ 92 


Shfh^ + 5A4( + /ift 

15 


'- 9 I 


—A 4 ^3/i^ + 3h‘^ + + 18 ^ 2 A 3 A 4 + A 4 + 3 A 5 + 3Ag + 3A7^ 

Ahfh'^ + — A4^17/i^ + 5/i^ + 15/i^^ + 12 A 4 ^Ai + A 2 + A 3 + 2 A 4 ^ 
+ 6(^8A2 + 7A2 + 4A6A7 + 7A2) 


9 

2^^ 

12/1? 


3hf + Sht + hi 


hlhl 


2A[hl + hl] + 24 A 3 + 33 A 4 


A 2 A 4 + 2 A 3 A 4 + A? + 2A? + A 6 A 7 + 5 A 7 


I2hl + Alt 


A 1 A 4 + 2 A 3 A 4 + A? + 2A? + 5A? + A 6 A 7 


7 A 4 

2Ai 

4 A 3 

2 A 4 


4 ( Ai + A2 


A? + A? 


37A? + IOA 6 A 7 + 5A 


IOA3A4 + 5 A? + 6 A? 


2 A 5 


5A? + IOA 6 A 7 + 37A? 


7 A 3 A 4 + 7A? + 12A? + 18A? + 2 OA 6 A 7 + 18A^ 


13A? + 34A? + 8 OA 6 A 7 + 34A^ 


16Af 


5A? + 6 AgA 7 + 5 A 7 


(A25) 


38 



As: 

/3^5^ = 2A5 Ai + A2 + 4A3 + 6A4 + lOAg + 4A6A7 + IOA7 — 2A5 7 i ^ + 72 ^ — 2A5 7^^^ + 72^^ 

(A 26 ) 

^2) , r 231 4 57 2 2 1413 4I 2^ r,2 ,2l 

/^As — As- —g2 + ^ 929 i + ~200 ^^ 405f3As + h^, 

4" ^92 ^As ^3/ij + 3h^ + h^^'j + 48As ^A 3 + 2 A 4 ^ + 72 ^Ag + A^^ 

+ ^As ^17^t 4“ 5/i^ + 15/i^^ — 24As ^Ai + A 2 — 8 A 3 — 12 A 4 ^ + 48 ^5Ag — AeAy + SA^^ 

— “As 3hf + 3h^ 4“ — hf 33h^\^ + 12 A 5 ^A 2 + 2 A 3 + 3A4^ + 12 A 7 ^A 6 + 5A7^ 

— ^ 12 /i^ + 4/i^^ As ^Ai + 2 A 3 + 3A4^ + Ag ^5Ag + A 7 ^ 

— 7As A 4 + A 2 — 4As Ai + A 2 IOA 3 + IIA 4 

- 2Ai [37A2 + 10AgA7 + 5A^j - 2 A 2 k\l + 10AgA7 + 37\‘^j 

- 4A3 [7A3A5 + I9A4AS + 18 A 2 + 20AgA7 + 18 A 2 

- 4 A 4 [ 8 A 4 A 5 + 19A2 + 22 A 6 A 7 + 19A2j + 6\l - 8 A 5 kxl + 2 IA 6 A 7 + 9A^ 

(A 27 ) 


39 



Ae: 

= 2 \q 6Ai + 3A3 + 4A4 + 5A5 + 2A7 3A3 + 2A4 + A5 — Ag 37I ^ + 72 ^ — Ag 37!^^ + 72 ^^ 

(A 28 ) 

~ 141 Ag + 9OA7 + —gldi 29 Ag + IOA7 + 187 Ag + 3OA7 + 205f3Ag + 3 hl 

+ —5^2 3 Ag + 3 /i^ + + 48 Ag ^Ai + Ag^ + 16 Ag ^Ag + 2 A 4 ^ + I6A7 ^2A3 + A4^ 

+ ^Ae + 15 /ig + 45 /i^^ + 48 Ag ^ 9 Ai + 3A3 + 5A4 + lOAg^ + 48A7 ^ 4 A 4 — Ag^ 

- ^Ag 27 ht + Sih'ihl + 33 ht + ll/i^ 

- Qlif Ag ^ 3 A 3 + 4A4 + 5 Ag^ + A7 ^6A3 + 4A4 + 2 Ag^ 

- ^6/ig + 2/i^^Ag 12Ai + 3A3 + 4A4 + 5 Ag 

- ^Ag 53 A? - A^ - 2 AiAg 33A3 + 35A4 + 37 Ag - 2A2A6 9A3 + 7A4 + 5 Ag 

- 2 Ag [l 6 A^ + 34A3A4 + 36A3A5 + 17 A^ + 38 A 4 Ag + 18 A^ 

- 2 (^Ai + A2) A7 9A3 + 7A4 + 5 Ag - 2A7 18 A^ + 28A3A4 + 20A3Ag + 17 A^ + 22A4Ag + 21A^ 

- 3 [ 37 A^ + 42 A^A 7 + llAgA? + 14 A? 

(A 29 ) 


40 



— 2A7 6A2 + 3A3 + 4A4 + 5A5 + 2 \q 3A3 + 2A4 + A5 — A7 7I ^ + 372 ^ 7^^^ + 372^^ 

(A30) 

= -^92 QOAe — 14 IA 7 + —g^gl lOAg + 29 A 7 + 30A6 + I 87 A 7 + 205 f 3 A 6 Sh^ + hi 

+ g 92 ^A 7 {^hl + 3 /i^ + + I44A7^A 2 + As^ + 48A3 ^ 2 Ag + A 7 ^ + 48A4^Ag + 2 A 7 ^ 

+ h\i(^51hl + hhl + 15hl^ + 48Ag ^ 6 A 3 + 4 A 4 — Ag^ + 48A7^9A2 + 3 A 3 + 5 A 4 + lOAg^ 

- -A7 llh^ + 28 hlhl + 9 hl + 3 /i^ 

- 6/1^2 A7 [12A2 + 3A3 + 4A4 + SAg 

- ^ 6 /ig + 2hl^ Ag ^ 6 A 3 + 4 A 4 + 2Ag^ + A 7 ^3A3 + 4 A 4 + 5Ag^ 

- ^A 7 -\l + 53A^ - 2 A 1 A 7 9 A 3 + 7 A 4 + 5Ag - 2 A 2 A 7 33 A 3 + 35 A 4 + 37 A 5 

- 2 A 7 16A^ + 34 A 3 A 4 + 36A3Ag + 17A^ + 38A4Ag + 18A^ 

- 2 (^Ai + A 2 ) Ag 9 A 3 + 7 A 4 + 5Ag - 2Ag 18A^ + 28 A 3 A 4 + 20 A 3 Ag + 17A^ + 22 A 4 Ag + 21 A^ 

- 3 14A^ + 11A^A7 + 42AgA? + 37A? 

(A31) 

[1] S. Chatrchyan et al. [CMS Collaboration], Phys. Lett. B 716, 30 (2012) [arXiv; 1207.7235 
[hep-ex]]. 

[2] G. Aad et al. [ATLAS Collaboration], Phys. Lett. B 716, 1 (2012) [arXiv:1207.7214 [hep-ex]]. 

[3] G. Aad et al. [ATLAS Collaboration], Phys. Rev. D 90, no. 5, 052004 (2014) [arXiv; 1406.3827 
[hep-ex]]. 

[4] V. Khachatryan et al. [CMS Collaboration], Eur. Phys. J. C 75, no. 5, 212 (2015) 

[arXiv: 1412.8662 [hep-ex]]. 

[5] G. Aad et al. [ATLAS and CMS Collaborations], Phys. Rev. Lett. 114, 191803 (2015) 

[arXiv: 1503.07589 [hep-ex]]. 

[ 6 ] S. Heinemeyer, O. Stal and G. Weiglein, Phys. Lett. B 710, 201 (2012) [arXiv:1112.3026 [hep- 
ph]], L. J. Hall, D. Pinner and J. T. Ruderman, JHEP 1204, 131 (2012) [arXiv:1112.2703 [hep- 


41 



ph]], U. Ellwanger, JHEP 1203, 044 (2012) [arXiv:1112.3548 [hep-ph]], P. Draper, P. Meade, 
M. Reece and D. Shih, Phys. Rev. D 85, 095007 (2012) [arXiv:1112.3068 [hep-ph]], A. Arbey, 
M. Battaglia, A. Djouadi, E. Mahmoudi and J. Quevillon, Phys. Lett. B 708, 162 (2012) 
[arXiv:1112.3028 [hep-ph]], M. Carena, S. Gori, N. R. Shah and C. E. M. Wagner, JHEP 
1203, 014 (2012) [arXiv;1112.3336 [hep-ph]]. 

[7] H. E. Haber and R. Hempfling, Phys. Rev. Lett. 66, 1815 (1991); Y. Okada, M. Yamaguchi 
and T. Yanagida, Prog. Theor. Phys. 85, 1 (1991); J. R. Ellis, G. Ridolfi and F. Zwirner, 
Phys. Lett. B 262, 477 (1991). 

[8] S. Heinemeyer, W. Hollik and G. Weiglein, Phys. Rev. D 58, 091701 (1998) [hep-ph/9803277]; 
S. Heinemeyer, W. Hollik and G. Weiglein, Eur. Phys. J. C 9, 343 (1999) [hep-ph/9812472]; 
S. Heinemeyer, W. Hollik and G. Weiglein, Phys. Lett. B 440, 296 (1998) [hep-ph/9807423]. 

[9] S. P. Martin, Phys. Rev. D 71, 016012 (2005) [hep-ph/0405022]. 

[10] G. Degrassi, P. Slavich and F. Zwirner, Nucl. Phys. B 611, 403 (2001) [hep-ph/0105096]. 

[11] A. Brignole, G. Degrassi, P. Slavich and F. Zwirner, Nucl. Phys. B 643, 79 (2002) [hep- 
ph/0206101]. 

[12] A. Dedes, G. Degrassi and P. Slavich, Nucl. Phys. B 672, 144 (2003) [hep-ph/0305127]. 

[13] S. P. Martin, Phys. Rev. D 66, 096001 (2002) [hep-ph/0206136]. 

[14] J. R. Espinosa and R. -J. Zhang, JHEP 0003, 026 (2000) [hep-ph/9912236]. 

[15] J. R. Espinosa and R. -J. Zhang, Nucl. Phys. B 586, 3 (2000) [hep-ph/0003246]. 

[16] S. P. Martin, Phys. Rev. D 65, 116003 (2002) [hep-ph/0111209]. 

[17] S. P. Martin, Phys. Rev. D 75, 055005 (2007) [hep-ph/0701051]. 

[18] R. V. Harlander, P. Kant, L. Mihaila and M. Steinhauser, Phys. Rev. Lett. 100, 191602 (2008) 
[Phys. Rev. Lett. 101, 039901 (2008)] [arXiv;0803.0672 [hep-ph]]. 

[19] P. Kant, R. V. Harlander, L. Mihaila and M. Steinhauser, JHEP 1008, 104 (2010) 
[arXiv: 1005.5709 [hep-ph]]. 

[20] J. L. Feng, P. Kant, S. Profumo and D. Sanford, Phys. Rev. Lett. Ill, 131802 (2013) 
[arXiv: 1306.2318 [hep-ph]]. 

[21] J. A. Gasas, J. R. Espinosa, M. Quiros and A. Riotto, Nucl. Phys. B 436, 3 (1995) [Erratum- 
ibid. B 439, 466 (1995)] [hep-ph/9407389]. 

[22] M. S. Carena, J. R. Espinosa, M. Quiros and C. E. M. Wagner, Phys. Lett. B 355, 209 (1995) 
[hep-ph/9504316]. 


42 



[23] M. S. Carena, M. Quiros and C. E. M. Wagner, Nucl. Phys. B 461, 407 (1996) [hep- 
ph/9508343]. 

[24] H. E. Haber, R. Hempfling and A. H. Hoang, Z. Phys. C 75, 539 (1997) [hep-ph/9609331]. 

[25] M. S. Carena, H. E. Haber, S. Heinemeyer, W. Hollik, C. E. M. Wagner and G. Weiglein, 
Nucl. Phys. B 580, 29 (2000) [hep-ph/0001002]. 

[26] P. Draper, G. Lee and C. E. M. Wagner, Phys. Rev. D 89, no. 5, 055023 (2014) 
[arXiv: 1312.5743 [hep-ph]]. 

[27] E. Bagnaschi, G. E. Giudice, P. Slavich and A. Strumia, JHEP 1409, 092 (2014) 
[arXiv: 1407.4081 [hep-ph]]. 

[28] J. P. Vega and G. Villadoro, JHEP 1507, 159 (2015) [arXiv:1504.05200 [hep-ph]]. 

[29] K. Cheung, R. Huo, J. S. Lee and Y. L. Sming Tsai, JHEP 1504, 151 (2015) [arXiv;1411.7329 
[hep-ph]]. 

[30] H. E. Haber and R. Hempfling, Phys. Rev. D 48, 4280 (1993) [hep-ph/9307201]. 

[31] J. E. Gunion and H. E. Haber, Phys. Rev. D 67, 075019 (2003) [hep-ph/0207010]. 

[32] J. Guasch, W. Hollik and S. Penaranda, Phys. Lett. B 515, 367 (2001) [hep-ph/0106027]. 

[33] M. S. Carena, J. R. Ellis, A. Pilaftsis and C. E. M. Wagner, Nucl. Phys. B 586, 92 (2000) 
[hep-ph/0003180]. 

[34] M. Carena, H. E. Haber, 1. Low, N. R. Shah and C. E. M. Wagner, Phys. Rev. D 91, no. 3, 
035003 (2015) [arXiv: 1410.4969 [hep-ph]]. 

[35] T. Ibrahim and P. Nath, Phys. Rev. D 63, 035009 (2001) [hep-ph/0008237]. 

[36] B. Li and C. E. M. Wagner, Phys. Rev. D 91, 095019 (2015) [arXiv:1502.02210 [hep-ph]]. 

[37] D. Buttazzo, G. Degrassi, P. P. Giardino, G. E. Giudice, E. Sala, A. Salvio and A. Strumia, 
arXiv: 1307.3536 [hep-ph]. 

[38] H. Arason, D. J. Gastano, B. Keszthelyi, S. Mikaelian, E. J. Piard, P. Ramond and 
B. D. Wright, Phys. Rev. D 46, 3945 (1992). 

[39] M. E. Machacek and M. T. Vaughn, Nucl. Phys. B 222, 83 (1983). 

[40] M. E. Machacek and M. T. Vaughn, Nucl. Phys. B 236, 221 (1984). 

[41] M. E. Machacek and M. T. Vaughn, Nucl. Phys. B 249, 70 (1985). 

[42] M. -X. Luo and Y. Xiao, Phys. Rev. Lett. 90, 011601 (2003) [hep-ph/0207271]. 

[43] L. N. Mihaila, J. Salomon and M. Steinhauser, Phys. Rev. Lett. 108, 151602 (2012) 
[arXiv: 1201.5868 [hep-ph]]. 


43 



[44] L. N. Mihaila, J. Salomon and M. Steinhauser, Phys. Rev. D 86, 096008 (2012) 
[arXiv: 1208.3357 [hep-ph]]. 

[45] K. G. Chetyrkin and M. F. Zoller, JHEP 1206, 033 (2012) [arXiv:1205.2892 [hep-ph]]. 

[46] A. V. Bednyakov, A. F. Pikelner and V. N. Velizhanin, Phys. Lett. B 722, 336 (2013) 
[arXiv: 1212.6829 [hep-ph]]. 

[47] K. G. Chetyrkin and M. F. Zoller, JHEP 1304, 091 (2013) [arXiv:1303.2890 [hep-ph]]. 

[48] A. V. Bednyakov, A. E. Pikelner and V. N. Velizhanin, Nucl. Phys. B 875, 552 (2013) 
[arXiv: 1303.4364 [hep-ph]]. 

[49] G. E. Giudice and A. Strumia, Nucl. Phys. B 858, 63 (2012) [arXiv: 1108.6077 [hep-ph]]. 

[50] P. S. B. Dev and A. Pilaftsis, JHEP 1412, 024 (2014) [arXiv:1408.3405 [hep-ph]]. 

[51] J. Beringer et al. [Particle Data Group Collaboration], “Review of Particle Physics (RPP),” 
Phys. Rev. D 86, 010001 (2012). 

[52] L. Maiani, A. D. Polosa and V. Riquer, Phys. Lett. B 724, 274 (2013) [arXiv:1305.2172 [hep- 
ph]]; A. Djouadi, L. Maiani, G. Moreau, A. Polosa, J. Quevillon and V. Riquer, Eur. Phys. J. 
C 73, 2650 (2013) [arXiv:1307.5205 [hep-ph]]; A. Djouadi, L. Maiani, A. Polosa, J. Quevillon 
and V. Riquer, arXiv:1502.05653 [hep-ph]. 

[53] M. Carena, 1. Low, N. R. Shah and C. E. M. Wagner, JHEP 1404, 015 (2014) [arXiv: 1310.2248 
[hep-ph]]. 

[54] G. Degrassi, S. Heinemeyer, W. Hollik, P. Slavich and G. Weiglein, Eur. Phys. J. C 28, 
133 (2003) [hep-ph/0212020]. M. Frank, T. Hahn, S. Heinemeyer, W. Hollik, H. Rzehak and 
G. Weiglein, JHEP 0702, 047 (2007) [hep-ph/0611326]. 

[55] T. Hahn, S. Heinemeyer, W. Hollik, H. Rzehak and G. Weiglein, Phys. Rev. Lett. 112, no. 
14, 141801 (2014) [arXiv: 1312.4937 [hep-ph]]; 

[56] E. Bagnaschi, et al. and Frensch, Felix and Heinemeyer, LHCHXSWG-2015-002, https: // 
cds.cern.ch/record/2039911. 

[57] M. D. Goodsell, K. Nickel and F. Staub, Eur. Phys. J. C 75, no. 1, 32 (2015) [arXiv: 1411.0675 
[hep-ph]]. 


44 



